Optical absorption in quantum dots: 
Coupling to longitudinal optical phonons treated exactly 
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Optical transitions in a semiconductor quantum dot are theoretically investigated, with emphasis 
on the coupling to longitudinal optical phonons, and including excitonic effects. When limiting to 
a finite number of m electron and n hole levels in the dot, the model can be solved exactly within 
numerical accuracy. Crucial for this to work is the absence of dispersion of the phonons. A suitable 
orthogonalization procedure leaves only m(m + l)/2 + n(n+ 1)/2 — 2 phonon modes to be coupled to 
the electronic system. We calculate the linear optical polarization following a delta pulse excitation, 
and by a subsequent Fourier transformation the resulting optical absorption. This strict result is 
compared with a frequently used approximation modeling the absorption as a convolution between 
spectral functions of electron and hole, which tends to overestimate the effect of the phonon coupling. 
Numerical results are given for two electron and three hole states in a quantum dot made from the 
polar material CdSe. Parameter values are chosen such that a quantum dot with a resonant sublevel 
distance can be compared with a nonresonant one. 
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I. INTRODUCTION 

Quantum dots (QDs) based on semiconductor struc- 
tures have received great attention for almost two 
decades now [1]. Whereas early experiments were only 
able to measure ensemble averages over many quantum 
dots, it is nowadays possible to address single quantum 
dots individually [2]. Photoluminescence and more re- 
cently absorption [3] has been measured with high spa- 
tial resolution, showing distinct lines related to individual 
dots (due to well width fluctuations in a quantum well). 
Using magneto-photoluminescence, detailed information 
on energy levels and phonon coupling could be extracted 
for single QDs made from II- VI semiconductor material 
[4]. This and related experimental work has enormously 
stimulated investigations on the interaction between the 
carriers on confinement levels and the surrounding po- 
larizable medium on a microscopic level. A proper un- 
derstanding of the related dephasing mechanisms might 
become important in view of future quantum computa- 
tion applications based on semiconductor quantum dots 
[5]. 

Carrier relaxation and dephasing in quantum dots cru- 
cially depends on the interaction with lattice vibrations 
(phonons). At first sight, scattering with longitudinal- 
optical (LO) phonons seems to be possible only if a sub- 
level spacing matches the LO phonon energy hujQ. In 
other (non-resonant) cases, scattering is expected to be 
impossible (or at least strongly reduced) . This is the so- 
called phonon-bottleneck problem which has been dis- 
cussed intensively in the literature [6]. The argument 
relies on the application of Fermi's golden rule which de- 
mands energy conservation in each individual process. 
However, electron-phonon interaction in quantum con- 
fined systems [7, 8] gives rise to a far more complicated 



picture including: Non-Markovian scattering, formation 
of electron-phonon bound states with level repulsion, and 
phonon satellites spaced at multiples of Tiloq. 

The coupling of quantum dot levels to acoustic 
phonons is responsible for another set of features. Due 
to their dispersion, the possible energy transfer covers a 
range from zero to a maximum value (typically, a few 
meV) which is related to the dot size. Consequently, 
phonon satellites appear here as broad bands surround- 
ing the zero-phonon line [9, 10]. The broadening of this 
line itself is a subject of intense research. Similar to 
the phonon bottleneck problem mentioned above, an ad- 
vanced theory [11] predicts a finite broadening even in 
cases where the next confined level is far away in energy 
compared to the maximal possible energy transfer. 

In the present work, we want to focus on the interac- 
tion with LO phonons. This is relevant for quantum dots 
made from polar material like CdSe, where the Frohlich 
interaction with LO phonons is much stronger than the 
coupling to acoustic phonons. Krai et al. [12] have 
dealt with the LO phonon bottleneck problem and found 
within a standard self energy approach (second order in 
the interaction) an appreciable broadening of the levels 
which persists even outside the exact resonance. It has 
further been argued that relaxation properties can be ob- 
tained from a convolution of the spectral functions of the 
discrete energy levels. 

Aiming at a non-perturbative treatment, a big advan- 
tage is the dispersionless nature of the LO phonons which 
allows a numerically exact treatment. We have derived in 
our previous paper Ref. 13 an efficient scheme to calculate 
eigenenergies and eigenvectors, which rests upon a uni- 
tary transformation of the phonons into a new set where 
only a few modes are coupled to the electronic degrees of 
freedom. We have calculated the electron spectral func- 
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tion and shown that it consists of a series of discrete 
delta functions. They are distributed around the bare 
level energies and at multiples of the LO energy (phonon 
satellites). Under resonance conditions, the eigenener- 
gies are still split in form of avoided level crossing. This 
is reminiscent to the Rabi splitting in a two-level sys- 
tem coupled to monochromatic photons (which formally 
replace the phonons). The self consistent second Born 
approximation for the self energy gives some gross fea- 
tures of the spectrum but has broadened levels instead 
of the closely spaced discrete lines in the exact calcula- 
tion. Thus, we had concluded that this approximation 
fails if phonons are coupled to discrete electronic lev- 
els. Recently, our numerically exact treatment was used 
to assess results obtained via the Davydov's canonical 
transformation [14]. 

Looking at the influence of phonon interaction on the 
interband transitions, one has to refer to a seminal paper 
by Schmitt-Rink and coworkers [7]. They have shown 
that the phonon interaction in general will increase with 
confinement (i.e. with reducing the dot size). For the 
Frohlich coupling to LO phonons, however, things are 
more subtle since here the difference in charge distri- 
bution between initial (valence band sublevel) and fi- 
nal state (conduction band sublevel) enters. Now, un- 
der strong confinement, the sublevel wave functions are 
getting more similar. Therefore, the optical transition is 
accompanied with practical no change in local charge dis- 
tribution, and the net LO phonon coupling is drastically 
reduced. 

In Rcf. [12], it has been claimed that the interband 
absorption spectrum can be taken as a convolution of 
the one-particle spectral functions for the electron and 
the hole. In the present work which extends Ref. [13] to 
a two-band many-level situation, and allows to calculate 
the linear optical response exactly, we will demonstrate 
the failure of this convolution approach. Indeed, it is 
missing the correct charge redistribution. 

Absorption in QDs was first discussed in Ref. [15] 
with emphasis on the non-adiabatic electron-phonon cou- 
pling, labeled as phonon-assisted transitions in Ref. [13]. 
In Ref. [16], a detailed comparison between features in 
QDs and phenomena from quantum optics was presented, 
again focusing on the electron-phonon interaction. Re- 
cently, the absorption spectrum of individual QDs has 
been discussed for intra-band [17] as well as for inter- 
band transitions [18]. In the latter work, the valence 
band levels were treated without phonon interaction - 
which leads to an even simplified form of the convolution 
approach. However, the resonance condition (matching a 
sublevel distance by Hloq) can be much easier obtained in 
the valence band since the sublevel spacing is here smaller 
than in the conduction band. 

The paper is organized as follows. In section II, the 
model and the Hamiltonian are presented, which includes 
both carrier-phonon and carrier-carrier (Coulomb) inter- 
action. In section III, we define and perform the uni- 
tary transformation which reduces the number of Bosonic 



modes coupled to the Fermionic states. This works if the 
phonons have no dispersion, and makes a numerical di- 
agonalization of the Hamiltonian feasible. Equations for 
the linear polarization in terms of exact eigenstates are 
derived in section IV, and the approximate convolution 
of electron and hole spectral function is given as well. In 
section V, numerical results for two prototype quantum 
dots made from the polar material CdSe are presented 
and discussed. Several details on the transformations and 
a list of material parameters are given in the Appendix. 



II. THE MODEL 

We start with the standard Hamiltonian which cou- 
ples the band states in a semiconductor to the lattice 
displacement. For the quantum dot, the confined elec- 
tronic states in the conduction (valence) band are given 
by Fermionic creation operators cj (uj). The lattice vi- 
brations are taken as longitudinal optical phonons with- 
out dispersion, hui , and represented by Bosonic opera- 
tors 6 q : 

h = £ M>&£&q + E {fi c h + e >h) (!) 
+ E(V fet -q) { M ^ c h + KA v ) ■ 

The coupling matrix elements (a = c, v) stem from 
the Frohlich interaction applied to the confinement wave 
functions (see Appendix C for details). Note that there 
is no phonon coupling between valence and conduction 
band states since the energy gap is far greater than the 
LO phonon energy, while the distance between confined 
levels in one band e" — ef (a = c, v) may be well in the 
range of the phonon energy. Electron spin is not included 
here since spin relaxation is a slow process compared to 
the spin-conserving electron-phonon interaction. Simi- 
larly, phonon-assisted transitions into the continuum of 
wetting layer states are not considered. For simplicity, 
we restrict ourselves to the interaction with bulk phonon 
modes, leaving the more precise picture of confined and 
interface modes for future investigations [20]. 

For treating optical transitions between the filled va- 
lence band and the empty conduction band, it is con- 
venient to switch from the conduction-valence-band de- 
scription used in Eq. (1) to the electron-hole picture. This 
is accomplished by the replacement 

cj^ej, v 3 ^h). (2) 
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Using v^-vi = Sjj — h\hj the Hamiltonian is rewritten as 
q j 

+ E(V b -*) { M ^ e h - KiA h ) ( 4 ) 
+ £3 + 

Note that the phonon interaction now carries a negative 
sign for the hole states compared to the electron states, 
which can be traced back to the different charge sign of 
the excitations. In Appendix A we show that the last 
line can be dropped by renormalizing the ground state 
energy. 

As we are focusing on the linear response, the Coulomb 
interaction leads to the formation of a single exciton only. 
The relevant interaction term is 

H c = - Y^ v m e k h \ h i e i i ( 5 ) 

klij 

and has to be added to Eq. (3). VkUj denotes the Coulomb 
matrix element between sublevel states (for details, sec 
appendix C). Under strong confinement conditions con- 
sidered here, the confinement states are only reshaped 
marginally [19], and the Coulomb interaction leads to al- 
most rigid shifts of the transition energies (exciton bind- 
ing energies). However, non-diagonal Coulomb couplings 
have to be included since they can be of similar order as 
the phonon couplings, and modify the oscillator strengths 
as well. 



III. REDUCING THE PHONON SUBSPACE 

The phonon coupling in Eq. (3) involves only certain 
combinations of phonon operators, e.g. J2q^^ji^i- ^ 
we consider a finite number (say to) of electron sublevels 
in the dot, and n hole sublevels, there are to 2 + n 2 such 
combinations (as argued above, phonon-assisted transi- 
tions between conduction and valence bands are absent). 
Since we are dealing exclusively with confined states, the 
confinement wave functions can be chosen to be real, and 
the symmetry ■ = holds. This reduces the lin- 

ear independent combinations to m(m+l)/2+n(n+l)/2. 

A further reduction can be achieved since in the 
present Hamiltonian Eq. (3), the number of electrons 
(N e ) and of holes (Nh) are conserved quantities. There- 
fore, two (diagonal) Fcrmion pairs in the interaction can 
be expressed by the remaining ones. We choose 

m— 1 n— 1 

eL e m = N e - e ] e 3 > h) n h n = N h -^h] hj , (6) 

3=1 ' 3=1 

which gives additions to all the other diagonal coupling 
terms, and a c number remainder. In order to shorten 



the subsequent writing, we introduce pair indices A = 
1, . . . , (N— 2) which combine either two electron sublevels 
(e:jl) or two hole sublevels (h:j7), with j < I. Caring for 
the different signs in the interaction of Eq. (3), we define 
in the conduction band 

= - Sj,M^ mm , (7) 

and in the valence band 

M qA = -M^.+^M^ n . (8) 

For treating the c number part properly, we have to in- 
troduce one further index A = k = N—l, with the matrix 
element 

Mq K = N e Mq mm — NhMq nn . (9) 
Now, the interaction term of Eq. (3) reads 

K-l 

J2(A X + A{) c) x + (A K + Al) (10) 

with a shorthand writing for the Fcrmionic operators. 
Note that for the nondiagonal terms j ^ I, we have to 
set e.g. 

(c t c) A = e]e J + e | ej . (11) 

The combinations of phonon operators entering Eq. (10) 
are 

Ax = J2 M <^ b i (A = l...«), (12) 
q 

but the Ax do not form an orthonormal set. As detailed 
in Appendix B, we apply the Gram-Schmidt orthonor- 
malization scheme to generate new phonon operators B\ 
which arc properly orthogonalized. The transformation 
is a linear one, 

A 

Ax = h a B a , (13) 

a=l 

and leads to the transformed Hamiltonian 

H = J2 ^oB{ Bx + Y (4 4 e 3 3 h) 

A 3 

+ EE(^aS a + /; Q st)( c t c ) A (14) 

A=l a=l 

K 

+ E (^B a + ILBl) . 

a=l 

To arrive at the standard diagonal form of the free 
phonon part, it was essential that the phonons have no 
dispersion. Otherwise, the subset of B x would still mix 
with the remaining phonon operators. These remaining 
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(a) (b) (c) (d) 

FIG. 1: Diagrammatic expansion of the optical polarization 
up to first order in the electron-phonon interaction. Lines 
denote the electron (e) and hole (h) propagator, the dotted 
line the phonon propagator, and the large dot stands for the 
optical dipole matrix element. The zeroth-order diagram (a) 
describes interband transitions between confined states un- 
affected by the phonon interaction. In the convolution ap- 
proach, only self-energy type diagrams (b, c) are considered, 
while in the full evaluation, vertex type diagrams (d) are in- 
cluded as well. 



degrees of freedom only contribute to the free phonon 
energy and are omitted from Eq. (14). 

A careful inspection of Eq. (14) shows that the Bosonic 
operator B K only appears as last element in the last line 
of Eq. (14), and in the free phonon part. It is therefore 
decoupled from the Hamiltonian and allows an indepen- 
dent solution in terms of a shifted oscillator, 

B K = B K + I KK /fuv ■ (15) 

Therefore, only 

m(m + l) n(n+l)_ 

2 2 K ' 

new modes are coupled to the Fermionic Hilbert space. 

Assuming a symmetric dot shape, the confinement 
functions are either even or odd, and the matrix elements 
M^ji have a definite parity in q, too. Consequently, the 
compound matrix elements M\ a defined in Eq. (B5) van- 
ish if A and a refer to different parity. Therefore, the 
tridiagonal matrix K (and I as well) has a block struc- 
ture, allowing the Gram-Schmidt procedure to work in 
each block independently. Note that the Coulomb ma- 
trix elements have an equivalent parity, and therefore do 
not mix these blocks. For certain confinement potentials 
like the harmonic potential, the number of phonon modes 
can be reduced even further (see appendix C). 

IV. LINEAR OPTICAL RESPONSE 

In this section, we want to contrast the direct evalu- 
ation of the time-dependent linear polarization and the 
absorption spectrum with the so-called convolution ap- 
proach: Here, only the spectral functions of the elec- 
tron and hole states are calculated, and the absorption 
is taken as a convolution of both quantities. Before de- 
riving the corresponding formal expressions, let us point 
out the main difference: For the direct evaluation, we 
start from the electron-hole vacuum (N e — 0, Nh = 0) as 



initial state and end up in the subspace of one electron- 
hole pair (N e — 1, Nh = 1). For the convolution ap- 
proach, quite different subspaces are invoked, namely 
N e = 0, N h = -> N e = 1, N h = for the electron spec- 
tral function, and N e = 0, N h = -> N e = 0, N h = 1 
for the hole one. In diagram language, any process which 
has a phonon correlation between electron and hole lev- 
els is discarded in the convolution approach. In Fig. 1, 
we display the relevant first order diagrams. The ver- 
tex diagram Fig. Id is the first one describing interband 
phonon correlations. In particular under strong confine- 
ment conditions, it compensates to a large extent the 
self-energy type diagrams (Fig. lb, c). Only the latter are 
kept in the convolution approach and we expect much to 
strong phonon satellites in this approximate treatment. 
Formally, in the convolution approach, the matrix ele- 
ments appear as (M c ) 2 or (M v ) 2 , while the full version 
includes the vertex corrections containing M C M V as well. 
For level-diagonal matrix elements, this combines into 
(M c — M v ) 2 . The appearance of matrix element differ- 
ences can be seen already in Eq. (3) . Therefore, the near 
cancellation of changes in the charge distribution [7] can 
only be achieved when treating self energy and vertex di- 
agrams on an equal footing. Note that excitonic effects 
are completely absent within the convolution approach. 

A. Direct evaluation 

The coupling to the light is described within the dipole 
approximation for interband transitions, adding 

E(t)Y^(jHj &j + h.c.) (17) 

ij 

as a perturbation to the Hamiltonian. The dipole matrix 
elements are given by an integral over the confinement 
functions, 

fMj =n cv jd 3 r^(r)^(r), (18) 

having as prefactor the dipole moment between the va- 
lence and conduction band, fx cv . 

The linear optical response follows from the polariza- 
tion after a unit-area delta pulse at t = and is given 
as dipole-dipole-correlation function. Expressed via the 
electron and hole operators, we have for t > 

P(t) = * E tHj{hi{t) e k (t) ej(0) fct( )) (19) 

with the time dependence of the operators in the Heisen- 
berg picture. The expectation value is shorthand writing 
for the statistical sum over initial states which contain no 
electron- hole pairs, and a thermal distribution of phonon 
mode occupations n\ — 0,1,.. . oo. This is denoted by 
\ri\, 0), with total energy 

E Q = huj E n x . (20) 

A 
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Then, we proceed with 
P(t) ^iJ2e- pE "J2^i^ 3 ( n ^0\h l (t)e k (t)\n x ,ij) . 

n x ij,kl 

(21) 

In practice, the independent summation over phonon oc- 
cupations n\ can be restricted to a maximum number 
which depends on coupling strength and temperature, 
f3 = 1/fegT. For simplicity, we have omitted a prefactor 
which ensures the proper normalization of the statistical 
sum. 

While the zero-pair states |nA,0) diagonalize the zero- 
pair Hamiltonian properly, we need to look for the one- 
pair states from the eigenvalue problem 

(H + H C ) 1 \4>)=E^). (22) 

Due to the reduction of the number of Bosonic modes, 
this can be solved numerically exact by expanding into 
the noninteracting one-pair basis \n\,ij). 

Plugging all time dependencies together, we obtain 

P{t) = iJ2^ E ° E Vki^j (23) 

"A ij,kl,if> 

x e* E ° ~ E ^' h {n Xl km(4>\n x ,ij). 
Starting with the initial value 

P(t = 0)^5>-^5> H | 2 , (24) 

nx kl 

the polarization evolves in time as a sum over many in- 
dividual oscillations. Although each of these terms does 
not have a damping, in the sum a general decay can be 
observed (quasi-dephasing) , as shown in Sec. V. 

The imaginary part of the Fourier transformation of 
the polarization function yields the absorption spectrum, 

/oo 
dtP{t)e lut (25) 
-OO 

« E e_/3Bo E ^iMj 

n x ij,kl,cf> 

x {n x ,kl\4>){il>\nx,ij)6(Et-Eo-hu>). (26) 

Let us point out that the exact absorption spectrum con- 
sists of delta peaks only. This was clear from the begin- 
ning since a finite perturbation cannot change the charac- 
ter of the spectrum of the unperturbed system. Since we 
started from dispersionless Bosonic modes, the discrete 
electronic spectrum cannot be altered by the non-zero 
electron-phonon interaction [21]. 
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FIG. 2: Schematic band diagram with sublevel spacings cor- 
responding to the resonant quantum dot. Vertical arrows - 
optical transitions, dotted lines - phonon processes. For fur- 
ther explanation, see text. 



Indeed, avoiding the evaluation of the two-particle (i.e. 
dipole-dipole) correlation functions would be an impor- 
tant reduction of the numerical labor, since the spectral 
function is a genuine one-particle function. We concen- 
trate on the electron spectral function A e ki (uj) which is 
defined by 

1 f°° 

A e ki (uj) = -Re dte^(e k (t)el(0)). (28) 
n Jo 

Here, we need the exact eigenstates of the Hilbert space 
with N e = 1, Nh = which will be denoted by H e \x) = 
E x \x)- Similar arguments as used in deriving Eq. (23) 
lead to 

A e u {u) ex £e"^° (29) 

x ^2(nx,k\x)(x\n\,i)S(E x - £ - M > 
x 

where \n\,k) is the (non- interacting) one-electron basis. 
Similar equations hold for the hole spectral function. Due 
to the assumed dot symmetry, the only off-diagonal spec- 
tral function will be the hole spectral function related to 
h:13. 



B. Convolution approach 

It has been argued in the literature [12], that the con- 
volution of electron and hole spectral functions may give 
a reasonable approximation to the absorption spectrum, 

a co »M oc Mh My AUu') Afou - w') . (27) 

ij,kl 



V. NUMERICAL RESULTS 

We consider a prototype quantum dot of ellipsoidal 
shape with CdSe as the semiconductor material. The 
relevant LO phonon energy is Tlujo — 24 meV. The larger 
mass of the holes results in a smaller level distance in 
the valence band. In order to cover a comparable en- 
ergy range in the conduction and the valence band, we 
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take into account the two lowest confined states in the 
conduction band (m — 2) and the three lowest confined 
states in the valence band (n = 3). In accordance with 
Eq. (16), we then have seven new phonon modes. 

According to what has been said before, there is the 
odd parity group (A = e:12, h:12, h:23) and the even 
parity group (A = e:ll, h:ll, h:22). The special mode A = 
k = 7 refers to the even parity group. Consequently, K 
and / decompose into a three-dimensional block for both 
parity groups. Due to the harmonic potential, no extra 
Bosonic mode is assigned to the transition h:13 which 
can be expressed by other Bosonic modes (see Appendix 
C). 

For characterizing the different dot sizes we define the 
first sublevel distance in each band (without phonon in- 
teraction and exciton effects) as 

A a ^\e a 2 -e1\ (a = c,v). (30) 




t[ps] 

FIG. 3: Femporal decay of the linear polarization for the res- 
onant QD at T = 77 K (upper curves) and T = 300 K (lower 
curves). The excitation pulse is centered around the lowest 
optical transition (hl-el exciton) with a spectral window of 
±12meV. The inset shows the initial decay. 



In the sequel, we calculate the linear optical properties 
for two selected dot sizes: The resonant quantum dot has 
perfect resonance between phonon energy and sublevel 
distance in the valence band, A" = Tiloq — 24meV, and 
consequently A c = 55meV. In Fig. 2, the confinement 
levels for the resonant QD are shown schematically. Ver- 
tical arrows mark the nonzero optical transitions. The 
dotted lines denote phonon processes. In the valence 
band, this is resonant with the sublevel distance (real 
transition, driven by M^ 12 ). In the conduction band, 
a virtual process involving the diagonal matrix element 
Mq U is depicted. 

The second choice is backed by experimentally mea- 
sured level distances from Refs. [4, 22] and comes out to 
be a non-resonant quantum dot with A 11 — 35 meV and 
A c = 80 mcV. More details on the extraction of size pa- 
rameters are given in Appendix D. 



A. Direct evaluation 

To determine the full linear response, we need to 
solve the eigenvalue problem for the one-pair subspace 
Eq. (22). As mentioned at the end of Sec. II, it will con- 
sist of six Fermionic pair states which are coupled to six 
Bosonic modes. We truncate the Hilbert space by al- 
lowing the sum of occupation numbers in all six Bosonic 
modes not to exceed seven. This gives 1716 different 
Bosonic states using the binomial coefficient (13 over 6). 
For the Fermionic states we can exploit an addi- 
tional parity symmetry. Thus, in total, two matrices with 
5148 x 5148 elements each have to be diagonalized. At 
T = 77 K, we have flhuio = 3.62, which tells us that LO 
phonon emission processes are dominant. By checking 
the convergence, we found that the first 1000 eigenvalues 
are sufficient to describe the spectra properly. 

In dealing with the special phonon mode -B7, we have 
to evaluate the matrix elements between the unshifted 
oscillator (in the initial state) and the shifted one (in the 



final state), 



(n,0\n',s)= (31) 

min(nV) 

e -s 2 /2 V n!n! „n+n'-2k 



k=Q 



k\ (n-k)\ {ri - k)\ 



The shift parameter is given as second argument in the 
state, here s = Ir?/hujQ. In the polarization, this part 
gives rise to an additional factor of 



nt,n' 7 



(32) 



Note the appearance of s 2 u>o as a correction to the final 
state energy. For the absorption spectrum, first the spec- 
trum without the special mode is calculated, and subse- 
quently spectrally displaced and added up, as the Fourier 
transform of Eq. (32) dictates. 

The temporal decay of the polarization amplitude 
\P(t)\ is displayed in Fig. 3 for the resonant QD. In 
the calculation, we have taken into account only tran- 
sitions which are energetically close to the (lowest) hl- 
el exciton. Thus, phonon satellites and the interference 
with the other transitions are suppressed. This spectral 
window would correspond to a finite-duration excitation 
pulse. For an elevated temperature of 300 K, the initial 
decay goes much faster, which resembles a traditional 
temperature-dependent dephasing. However, at larger 
times, the polarization oscillates irregularly around a fi- 
nite value (therefore, we would like to use the term quasi- 
dephasing). 

A similar dynamical decay due to incommensurate en- 
ergies has been discussed in Refs. [16, 23]. In real sys- 
tems, these LO phonon beats will be damped finally by 
anharmonicity effects [24]. 

For the absorption spectrum, we have chosen to 
broaden all discrete lines in Eq. (25) with a fixed Gaus- 
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FIG. 4: Normalized absorption spectrum a(u>) for the non- 
resonant QD (left) and the resonant QD (right) at T = 77 K. 
The lower panels refer to the full model while the upper panels 
are restricted to level-diagonal phonon coupling. The bare 
hi— >el transition energy is chosen as zero of energy. 
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FIG. 5: Normalized absorption spectrum a(ui) for the non- 
resonant QD at T = 77 K. Upper panel: Without Coulomb 
interaction. The transition energies without phonon inter- 
action are shown as vertical dashed lines which refer to the 
vertical arrows of Fig. 2. Middle panel: All couplings and 
levels included. Lower panel: Only two hole levels included. 



sian of variance a = 1 meV. In this way, we can visu- 
alize both, the transition energies and their oscillator 
strengths. 

In Fig. 4, the absorption spectrum a(u>) is shown for 
both the non-resonant QD (left) and the resonant one 
(right). In the lower panels, all phonon couplings are in- 
cluded. A set of closely spaced levels (bundle) is seen 
around the two main optically allowed exciton transi- 
tions hl-el and h2-e2. The excited-state bundle for the 
resonant QD exhibits a particular broad range. A sim- 
ple theory like e.g. the second Born approximation [12] 
would give here a real broadening of the line, which may 
be interpreted as a real dephasing process via energy- 
conserving phonon emission. We see that in the present 
exact treatment, the situation is a bit more complex, and 
can at best be called quasi-dephasing. Also, the compar- 
ison with the non-resonant dot (left) shows that quasi- 
dephasing is in no way restricted to exact resonance. Al- 
though a treatment of exciton occupation is outside the 
frame of the present paper, this finding points to the 
absence of a clear phonon bottleneck in polar quantum 
dots. 

We have also calculated the spectrum for the reduced 
model of level-diagonal phonon coupling M^- cx 5y (up- 
per panels in Fig. 4) . This case is known to be ex- 
actly solvable and called Independent Boson Model, see 
Ref. [25]. We found excellent agreement supporting our 
method, and could justify the chosen truncation in Bo- 
son occupation. The level-diagonal absorption spectrum 
exhibits phonon satellites which carry only a few percent 
of the total weight, leaving the zero phonon line as the 
dominant feature. The inclusion of the inter-level or non- 
adiabatic coupling thus changes the absorption spectrum 
significantly which was already emphasized in Ref. [15]. 

In Fig. 5 we show the different contributions separately. 



Without Coulomb interaction (upper panel), the transi- 
tions are slightly shifted compared to their bare posi- 
tions marked by dashed vertical lines. This allows to ex- 
tract a polaron shift of the lowest transition in the non- 
resonant QD of 4meV. This is much smaller than the 
standard expression for the polaron shift in bulk CdSe 
(electrons: 11.3 meV, holes: 21.1 meV) and shows indeed 
the strong charge cancellation in a QD as discussed in 
Ref. [7]. In addition, the restriction to a finite number 
of sublevels might underestimate our calculated QD po- 
laron shift. The comparison with the full calculation in- 
cluding phonon and Coulomb interaction (middle panel) 
shows the importance of the exciton effect. There is a 
general shift which is related to the dominant diagonal 
Coulomb matrix element fun = 67.4meV. However, the 
competition between nondiagonal Coulomb interaction 
and phonon coupling (which are of the same order of 
magnitude) modifies the spectrum in a complex manner, 
and cannot be reduced to a rigid spectral shift. This 
clearly shows the importance of excitonic effects in the 
absorption spectrum of strongly polar quantum dots. 

In the lower panel, results of a calculation with only 
two hole levels are shown. The h3-el transition and its 
excitonic modification is, of course, missing here. Never- 
theless, the basic features in the spectrum prevail. The 
polaron shift is a bit smaller and the spectrum is less 
"smeared out" compared to the full model with three 
hole levels (middle panel). 



B. Convolution approach 

Performing the unitary transformation described in 
section III on the Hamiltonian for each band separately, 
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FIG. 6: Normalized absorption spectrum a c ° n (ui) obtained 
from convoluting the spectral functions for electron and hole 
in the non- resonant QD. The transition energies without 
phonon interaction are shown as vertical dashed lines. The 
upper panel refers to the diagonal phonon coupling, while the 
lower panel includes all coupling matrix elements. 



we end up with a model containing three Bosonic modes 
for the conduction band and six Bosonic modes for the va- 
lence band. Using the conservation of the electron (hole) 
number and again a parabolic confinement potential, 
we finally have to diagonalize a Hamiltonian with two 
Bosonic modes and two Fermionic states for the conduc- 
tion band and four Bosonic modes and three Fermionic 
states for the valence band. For the two-state model, 
this has been discussed at length in our previous paper 
Ref. [13]. 

Within the convolution approach, excitonic effects can 
only be included approximately by a rigid shift of the 
spectrum using an ad-hoc exciton binding energy. In the 
following discussion, we want to demonstrate that even 
the electron-phonon interaction is not treated appropri- 
ately. Therefore, we have no Coulomb interaction in this 
subsection. 

In Figure 6, the absorption spectrum a con {ui) from the 
convolution of the spectral functions Eq. (27) is shown for 
the non-resonant case. In order to reach the same (artifi- 
cial) broadening of the spectrum, we have to broaden the 
spectral functions with a Gaussian of reduced variance, 
a = l/V2meV. 

The inclusion of nondiagonal phonon matrix elements 
(lower panel) is not so important here, as the spectrum 
is dominated by strong phonon satellites. A compar- 
ison with the exact spectrum (neglecting excitonic ef- 
fects) of Fig. 5 (upper panel) shows, however, that these 
satellites are much too strong. As discussed at the be- 
ginning of Sec. IV, it is the correlated phonon scattering 
between electron and hole states which reduces the satel- 
lite structure appreciably. Thus, the one-particle spectral 
functions are not able to describe the optical transitions 
properly. Thus, we conclude that the convolution ap- 
proach as applied in Refs.[12, 18] is insufficient, even if 



the phonon coupling is of moderate strength, as in the 
present example. 



VI. SUMMARY 

We have presented a solvable model which describes 
the optical properties of a single quantum dot inter- 
acting with LO phonons, and includes exciton forma- 
tion. The linear polarization and the corresponding ab- 
sorption spectrum is calculated by diagonalizing the ap- 
propriate Hilbert space of six electronic levels and six 
Bosonic modes split off from the full phonon modes via 
an orthogonalization procedure. The correct evaluation 
of the two-particle dipole-dipole correlation function is 
contrasted with a simplified approach using the spectral 
convolution of one-particle spectral functions. The nu- 
merics shows this convolution approach to be inadequate, 
as it gives phonon satellites strongly enhanced compared 
to the correct result. Qualitatively, this can be traced 
back to a missing compensation between different contri- 
butions in a diagrammatic analysis. 

Parameter values for CdSe quantum dots are used, and 
a comparison between resonant and non-resonant QDs 
is carried out. Signatures of the phonon bottleneck are 
found, i.e. a stronger quasi-dephasing for the resonant 
case, but this is not as dramatic as simple arguments us- 
ing Fermi's golden rule for phonon emission would pre- 
dict. 
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APPENDIX A: ELECTRON-HOLE 
TRANSFORMATION 

Switching from the conduction-valence-band descrip- 
tion to the electron-hole picture produces a shift in the 
Bosonic operators. Here we show that only the shift of 
the zero-momentum mode survives. 

To proceed we have to use the complete form of the 
confinement wave function including the Wannier basis 
w a (v) of the band under consideration, 



*?(r) 



R 



R) 



(Al) 



where ?/)? R is the envelope part (called confinement func- 
tion in the previous sections). The orthogonality of the 
(real) wave functions induces 



$n = 



d 3 i 



*?(r)*?(r)= 5>? R ^? R 



R 



(A2) 
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since the Wannier functions are orthogonal on different 
lattice sites R^R'. Dropping prcfactors, we evaluate 

EKn « / *r E ^R^R'^ qr (A3) 

j Jn RR'j 

x w a (r - R) w a (r - ~R!) = C q . 

The sum over j gives a Kronecker symbol <5rr' due to 
the completeness of the coefficients V'jR- The remaining 
integral simplifies to 

c * = I ^E e4qru; '( r - R ) ( A4 ) 

JO R 

= / riVe^^E 6 ^ 
Jn d 



R 



The last sum over R produces 7V<5 q o, where A/" is the 
number of elementary cells O within the normalization 
volume, £1 = TVSIq- Altogether we have 



C q = <5q0 A/" 



(A5) 



Therefore, the correction in the last line of Eq. (3) reduces 
to zero momentum, 

(bl + & ) y/Afhuo a v (q) , (A6) 

In the prefactor a v (q), we let q — > at the end. 

In order to remove the last line from the Hamiltonian, 
the zero-momentum phonon operator is shifted according 
to 



bo + a v {q)^/Af/huo a — > b , 



(A7) 



which brings from the free phonon Hamiltonian a 
quadratic contribution J\fa^(q). Together with the va- 
lence energy sum, this gives an unimportant c num- 
ber contribution which can be dropped. However, the 
electron-phonon interaction is getting an addition at 
q — > 0, too. Due to the orthogonality of the confinement 
functions, we have in leading order 



M a , 



-^a v (q)(5 jl+ 0(q)) , (A8) 



which gives as correction 



-2a 2 v (q) (eje, - h\hi) . 



(A9) 



However, this term vanishes for any state having no or an 
equal number of electrons and holes (charge neutrality). 
For optical excitation and recombination, this is just the 
relevant sector of the Hilbcrt space. Therefore, the sin- 
gularity in a v (q) oc 1/q (which would have to be treated 
carefully) does not contribute to physical quantities. 

Finally, we note that the above argument also holds 
for acoustic phonons. 



APPENDIX B: GRAM-SCHMIDT 
ORTHO-NORMALIZATION 

We apply the Gram-Schmidt scheme by first construct- 
ing the orthogonal operator set 

A-l 

Bi = A lt B X = A X -Y / [Ax,Bl] B a , (Bl) 



and normalizing it afterwards via 
B x = N X B X where N7 2 = 



(B2) 



The new operators obey the canonical commutation rules 
[B x ,Bl]=5 Xa . (B3) 
The final result can be written as a linear transformation 



B X =N X J2 K Xa A a . 



(B4) 



a=l 



By construction, K Xa is a tridiagonal matrix which has 
nonzero elements for A > a only. The key ingredient for 
its evaluation is the commutator 

[A x ,Al] = E M qAM*, Q [fe q ,6 q ; 

qq' 

= ^M qA M* Q = M Aa , (B5) 
q 

and we get a recursive determination according to 

K xx = 1, (B6) 

A-l 

K Xa = -^NpKp^M^K^ (A>a). 

[j=a 



The norm follows from 



N 



x 2 = E -K Xa M a fjK X p . 

a, /3=1 



(B7) 



In order to transform the Hamilton operator to the new 
phonon operators B x , we need to invert Eq. (B4), 



A * = E haBa ■ 



(B8) 



a=l 



Again, I is tridiagonal and can be determined quite easily 
via 

A-l 

I xx = 1/N X , \>a: /a« = - E K ^ ■ ( B9 ) 

0=a 
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APPENDIX C: COUPLING MATRIX ELEMENTS 

The standard Frohlich coupling for the electron-LO- 
phonon interaction is applied to the dot confinement 
states, 



M% = 



£S 



(CI) 



fl 2q 2 eo \e, 

= J d 3 r^(r)e^iP?(r) (a = c,v). 

In the subspacc of one electron-hole pair states, the 
Coulomb interaction H c reduces to 



(n x ,kl\H c \n' x ,ij) = -5 nx y v kUj 



(C2) 



which is diagonal in the phonon quantum numbers. The 
matrix element is given by 



Vklij = Jd 3 rd 3 r'r k (v)^(r') 



Ane ss\r - r 



^VJ(r')^ c (r) 
(C3) 



Diagonalizing the phonon-free part together with H c 
would lead to (non-polar) exciton transition energies. 
However, in the present context it is more appropriate 
(and easier) to diagonalize phonon and excitonic effects 
together. Consequently, the main interband transition 
energies are related to exciton-polaron states in the quan- 
tum dot. 

For simplicity we consider an anisotropic parabolic po- 
tential as dot confinement for both, electrons and holes, 
with x as the long axis. The three energetically lowest 
wave functions are given by 



V?(r) 



i 1 / x 2 y" z" 

"~2 [x 2 + yJ + Z 2 



N, 



■ exp 



VSto - V2— <(r), 



(C4) 



where X a > Y a , Z a are the spatial extensions (variances) 
of the ground state, and N 2 — n 3 / 2 X a Y a Z a its normal- 
ization. 

The matrix elements in Eq. (CI) and Eq. (C3) thus 
read 

Ku = CX P (- \ (ll X l + iX + «X 
$q21 = iqxXa/yfiQ*!!, 

<fq 22 = (l-q 2 x X 2 a /2)^ u , (C5) 
<fq3 2 = * (fe^a - qlXl/4) $« n , 
*q 33 = (1 - &Xl + Q$Xi/8) $« n . 



Notice that $q 22 



%3i/V2. Therefore, one 
additional Bosonic mode can be eliminated. 

The new coupling constants M\ a are obtained by 
integrating a pair of coupling constants M° tJ over q, 
Eq. (B5). This final integration can be reduced to the 
following (elliptic) integrals, 



dt 



r 



Vt(t-p)(t-q) 



(n = 0,1,2,3,4) (C6) 



with p = 1 - (Y 2 + Y 2 )/(X 2 + X 2 ) < 1 and q = 1 - 
{Zl + Zl)j(X 2 a + Xl) < 1. For a cylindrically symmetric 
potential with Y a = Z a , we have p = q, and the integrals 
reduce to simple analytic functions, 

1 . 1 + VP . ^-(n+l/2)- 1 

Jo = — log , Jn+l = 



Vp 1 - Vp 



P 



(C7) 

The same type of integration over q appears in the 
Coulomb matrix elements Eq. (C3). For example, the 
diagonal matrix element for the lowest transition hl-cl 
can be evaluated as 



Jn=0 



(C8) 



where p and q in Eq. (C6) are determined with a = c, b = 
v. 



APPENDIX D: MATERIAL PARAMETERS 

Parameter values for the polar semiconductor material 
CdSe forming the quantum dot are from Ref. [26] as LO- 
phonon energy hui = 24mcV, conduction band mass 
m c = 0.13 mo, and valence band mass m v — 0.45 mo- 
Dielectric constants have been taken from Ref. [27] : Ss = 
9.57, £oo = 6.27. 

We assume that the confinement potentials of electrons 
and holes are scaled by a fixed ratio p which equals the 
ratio between conduction band offset and valence band 
offset for barrier (ZnSc) and dot material (CdSe), p = 
5.44 [26]. Applying this to the parabolic potential, we 
find a unique ratio between the lengths in the oscillator 
eigenfunctions Eq. (C4), 



m c p 



1/4 



, (L = X,Y,Z) (Dl) 



which equals £ = 0.89. 

Since the x direction is taken as the longest one, the 
energetic distance between ground level (j = 1) and 
first excited level (j = 2) is exclusively given by the x- 
confinement, 

(a = c,v). (D2) 



m a X 2 



Taking experimental values which have been reported 
for a certain CdSe quantum dot in Refs. [4, 22] (A c = 
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80 meV, A" = 35 meV) we obtain according to Eq. (D2) 
X c = 2.71 nm and X v = 2.20 nm. This yields a length 
ratio of £ = X c /X v = 1.23 which is larger than the value 
derived from Eq. (Dl). However, a possible alloying of 
the QD material [20] and deviations from the (idealized) 
parabolic confinement potential may introduce substan- 
tial uncertainties. In the numerical calculations, we have 
used £ = 1.23 throughout. 

There is no direct experimental access to the size of 
the QD in the other (shorter) directions. To keep things 
simple, we take the same representative value for the two 
smaller lengths: Y v = Z v = 1.3 nm and Y c = Z c = £Y V = 
1.6 nm. The next levels corresponding to this shorter 
length would have a spacing of 100 me V in the valence 
band. Since this is larger than twice the level spacing 
corresponding to the larger X extension (70meV), we 
have properly selected the lowest hole confinement states 
in Eq. (C4). 

It is easily seen from the parity of the oscillator wave 
functions Eq. (C4), that the optical transitions hl^e2, 
h2^el and h3^c2 arc forbidden. The remaining non- 
zero dipole matrix elements Eq. (18) can be expressed by 
the same length ratio £ introduced above. In units of /t c „, 
we find 



M31 



2£ 



■e 



1/2 



M22 = 



2£ 



3/2 



1/2 



(D3) 



With £ = 1.23, we get /in = 0.99, /i 22 = 0.97 and 
/t3i = 0.14. Such values very close to a Kronecker 
delta are typical for the strong confinement in a quan- 
tum dot. The energetic distance between the dominant 
transitions equals A c + A" = 115 meV (without polaron 
and Coulomb corrections). 



Apart from the parameter set deduced from the ex- 
perimentally given level distances, we chose another one 
which has perfect resonance between the hole level dis- 
tance and the LO phonon energy, i.e. A v — 24meV. 
This yields a somewhat larger QD with X v = 2.66 nm 
and X c = £X V = 3.27 nm, while keeping all the other pa- 
rameters the same. The level distance for the electrons 
amounts to be A c = 55 meV, giving an energetic distance 
between the main allowed transitions of 79 meV. We call 
this the resonant QD, while the parameter set backed by 
experiment is referred to as the non-resonant QD. 



The diagonal Coulomb matrix elements are in the 
range of 50-70 meV which is four times the exciton bind- 
ing energy in bulk CdSe (15meV) and quantifies the 
strong enhancement of Coulomb effects in small QDs. 
The non-diagonal Coulomb matrix elements are around 
15meV which is comparable to the effective phonon cou- 
pling strengths. 
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